Gravitational lensing with f(x) = x 3 ^ 2 gravity in accordance with astrophysical 

observations 
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In this article we perform a second order perturbation analysis of the gravitational metric theory 
of gravity f(x) = X 3//2 developed by Bernal et al. We show that the theory is capable to account 
exactly for two observational facts: (1) the phenomenology of flattened rotation curves through the 
Tully-Fisher relation observed in spiral galaxies, and (2) the details of observations of gravitational 
lensing in galaxies and groups of galaxies, without the need of any dark matter. We show how all 
dynamical observations on flat rotation curves and gravitational lensing can be synthesised in terms 
of the empirically required metric coefficients of any metric theory of gravity. We construct the 
corresponding metric components for the theory presented at second order in perturbation, which 
are shown to be perfectly compatible with the empirically derived ones. It is also shown that, 
in order to obtain a complete full agreement with the observational results, a specific signature 
of Riemann's tensor has to be chosen. This signature corresponds to the one most widely used 
nowadays in relativity theory. Also, a computational program, the MEXICAS (Metric EXtended- 
gravity Incorporated through a Computer Algebraic System) code, developed for its usage in the 
Computer Algebraic System (CAS) Maxima for working out perturbations on a metric theory of 
gravity is presented and made publicly available. 

PACS numbers: 04.50.Kd,04.20.Fy,04.25.Nx,95.30.Sf,98.80. Jk,98.62.Dm,98.62.Sb 
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I. INTRODUCTION 

When Einstein introduced his theory of general rela- 
tivity, an astrophysical prediction for the motion of the 
planet Mercury (a massive particle) through its orbit was 
made [l3| . The second step was to test general relativity 
through the deflection of light (massless particles) coming 
from stars appearing near the Sun's limb during a solar 
eclipse [12]. Both observations constituted the first co- 
herent steps towards the solid foundation of general rela- 
tivity, a theory capable of describing gravitation through 
a correct relativistic description. 

In this sense, any metric theory of gravity must be 
compatible with both kinds of observations, the dynami- 
cal ones for massive particles and the observations of the 
deflection of light for massless particles. The correct ap- 
proach is extensively described in the monograph by Will 
[44| where it is shown that when working with the weak 
field limit of a relativistic theory of gravity in a static 
spherically spacetime, the dynamics of massive particles 
determine the functional form of the time component of 
the metric, while the deflection of light determines the 
form of the radial one [see alsol45l. and references therein]. 

To order of magnitude and through a first perturbation 
analysis, Bernal et al. [|[ have shown that it is possible 
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to recover flat rotation curves and the Tully-Fisher rela- 
tion (i.e. a MONDian-like weak field limit) from a metric 
theory of gravity, which includes the mass of the system 
in the gravitational field's action. Such limit is of high 
astrophysical relevance at the scales of galaxies, where 
MOND accurately describes the rotation curves of spiral 
galaxies and the Tully-Fisher relation without the need of 
dark matter [see e.g. 0, |2j|. In this article we show the 
strength of the calculations made by Bernal et al. [H by 
doing an extensive analysis from perturbation theory for 
a static spherically symmetric metric and show that in 
the weak field limit our results are in perfect agreement 
not only with the Tully-Fisher relation, but are also in 
exact accordance with observations of gravitational lens- 
ing at a wide range of astrophysical scales. 

Extensions to Einstein's general theory of relativity 
have been pro pos ed since the publication of the theory 
itself [see e.g. [36J. However, it has not been until re- 
cent times that observations at different mass and length 
scales have concluded that in order to keep Einstein's 
field equations valid, unknown dark matter and energy 
entities need to be added to the theory. In this article, 
an equally important approach is taken where the ex- 
istence of these dark unknown entities is not required. 
We show the theory built by Bernal et al. [H to be in 
accordance not only with the very well established ob- 
servations of the dynamics of massive particles through 
the Tully-Fisher relation, but also with the dynamics of 
massless particles through the bending of light as astro- 
physically observed. 

Mendoza and Rosas-Guevara [28[ and Rosas-Guevara 
1351 showed for the first time that metric theories of 
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gravity are capable of producing more deflection of light 
that the one produced by Einstein's general relativity. 
This was done using the metric theory of gravity con- 
structed by Sobouti (38[. The implications of this re- 
sult invalidated the so called no-go theorem for metric 
f(R) theories of gravity proposed by Soussa [39(, Soussa 
and Woodard [40(. Furthermore, in the present work we 
show that it is possible to explain the observed gravita- 
tional lensing for galaxies, and groups of galaxies with- 
out the need of invoking dark matter. Developments by 
Capozziello, Cardone, and Troisi [see e.g. J5fl, Horvath 
et al. [see e.g. |23 |. Nzioki et al. [see e.g. Hj] in weak 
and strong lensing regimes of extended metric theories of 
gravity have followed the work by Mendoza and Rosas- 
Guevara [HI but are not completely coherent with dif- 
ferent astrophysical observations. 

Testing any metric theory of gravity against observa- 
tions can be cumbersome. From an action principle one 
must derive field equations, which in principle, have to 
be solved for e.g. in spherically symmetric spacetimes. 
The solutions to this lead to metric coefficients which 
in turn yield orbits for massive and massless particles, 
to be then compared against astrophysical observations. 
These last are varied and diverse e.g., centrifugal equi- 
librium orbits at a variety of radii, for systems having 
total masses spanning several orders of magnitude, and 
the observed shears and caustic positions of gravitational 
lensing observations. 

Fortunately, we have derived a much more direct and 
generic approach. First, dynamical observations regard- 
ing the amplitudes of galactic flat rotation curves satisfy 
a well known scaling with the fourth root of the total 
baryonic content: the Tully-Fisher relation. To second 
order in perturbations of the velocity measured in units 
of the speed of light, this can be shown to imply a defini- 
tive empirical form for the time component of any met- 
ric theory not requiring dark matter. Second, we show 
that all gravitational lensing observations on elliptical 
and spiral galaxies, as well as for groups of galaxies can 
be synthesised as the requirement for the same isother- 
mal total matter distribution as needed to explain the 
observed spiral rotation curves and dynamics about ellip- 
tical galaxies, if one assumes Einstein's general relativity. 
This in turn, from studying directly the lens equation in 
general relativity, implies a bending angle which is inde- 
pendent of the impact parameter, and which scales with 
the square root of the total baryonic mass of a system. 
It can then be shown that this, in combination with the 
empirical time component mentioned above, leads to a 
definitive empirical form for the radial component, for 
any metric theory not requiring dark matter. 

Thus, we synthesise all dynamical and gravitational 
lensing astrophysical observations at galactic and galaxy 
group scales, to second order in perturbation, into em- 
pirical time and radial metric components. It is through 
comparing the above to perturbed metric coefficients to 
the same order coming from the metric theory treated in 
this paper that we are able to show its full compatibil- 



ity with all relevant dynamical and gravitational lensing 
astrophysical observations. 

The article is organised as follows. In section [TT| the 
concept of weak field limit for a static spherically sym- 
metric metric is established and we define the relevant 
orders of perturbation to be used throughout the article. 

In section HTT1 we perturbe the vacuum field equations 
for the metric theory built by Bernal et al. [3j and show 
that for a point mass source they closely resemble the 
ones usually adopted in f(R) gravity in vacuum. How- 
ever, these equations slightly differ under the approxima- 
tions of the mass and length scales associated to galaxies 
and groups of galaxies -where gravity is expected to differ 
from Einstein's general relativity in the absence of any 
dark matter component. In section HVl we obtain the so- 
lution for the Ricci scalar up to the second order from 
the perturbed field equations and discuss the importance 
of the signature in the Riemann tensor to yield the cor- 
rect results. In section [V] we obtain the coefficients of 
the metric up to the second order in perturbation. In 
section I VII we obtain the metric coefficients up to the 
second order in an empirical way, without reference to 
any specific metric theory of gravity, using the dynami- 
cal phenomenology of galaxies and groups of galaxies and 
the gravitational lensing produced by these objects. In 
that section we also compare the metric coefficients ob- 
tained in[V]with those empirically obtained and show full 
consistency. Finally in section[VlT]we discuss our results. 

II. THE WEAK FIELD LIMIT 

An excellent account of perturbation theory applied to 
metric theories of gravity (in particular general relativity) 
can be found in the monograph written by Capozziello 
and Stabile Q, Will [H|. In this section we define the 
relevant properties of the perturbation theory having 
in mind applications to the metric theory developed by 
Bernal et al. [3[. 

Let us consider a fixed point mass M at the centre of 
coordinates generating a gravitational field. Under these 
considerations, the spacetime is static and its spherically 
symmetric metric g^ v is generated by the interval 

ds 2 = g^dx^dx" = g 00 c 2 dt 2 + gil dr 2 - r 2 dQ 2 . (1) 

In the previous equation and in what follows, Ein- 
stein's summation convention over repeated indices is 
used. Greek indices take values 0,1,2,3 and Latin ones 
1, 2, 3. As such, in spherical coordinates (x°, x 1 , x 2 , x 3 ) = 
(ct, r, 9, if), where c is the speed of light, t is the time co- 
ordinate, r the radial one, and 9 and ip are the polar 
and azimuthal angles respectively. Also, the angular dis- 
placement cir2 2 := d6> 2 +sin 2 9 dip 2 . Due to the symmetry 
of the problem, the unknown functions goo and g\ i are 
functions of the radial coordinate r only. We choose a 
(+, — , — , — ) signature for the metric of the spacetime. 
The radial component of the geodesic equations 
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ds 2 



dx 13 dx x 
ds ds 



= 0, 



(2) 



for the metric ((T|) in the weak field limit, i.e. when the 
speed of light c — ¥ oo, is given by 



1 d 2 r 
c 2 di 2 



1 



500,; 



(3) 



where the subscript ( := <9/dr denotes the deriva- 
tive with respect to the radial coordinate. In the above 
relation we have assumed that for the weak field limit 
ds = c dt and since the velocity v <§C c then v l <§; dx° /dt 
with v l := (dr/dt, rdO/dt, rsmOdip/dt). 

In this limit, a particle bounded to a circular orbit 
about the mass M experiences a centrifugal radial accel- 
eration given by: 



d 2 r 
ch 2 " 



v 
r 



(4) 



for a circular or tangential velocity v. 

When material particles are used as test particles in 
the weakest limit of the theory, the metric takes the form 
[see e.g. Hl| : 



ff 00 =ff 00(0) + 5 00(2)^ 1 _ g (2) +o(4); 
11(2) 



9 11 =.9 11(0) 

g 22 =5 22(0) 
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1-^+0(4), 



-1/r 2 , 



(7) 



,33 



=g 22 /sin 2 0. 



In fact, to the lowest order of perturbation, we need to 

(2) • (2) 

find the time <7q and radial g xl metric components up 
to the second order to compare with the astrophysical 
observations of material particles and bending of light 
MM- 



III. FIELD EQUATIONS 

Bcrnal et al. 3] proposed an extended gravitational 
field's action in the metric approach given by: 



16nGL 2 M 



/(x)v /3 5 d4 a ; 5 



(8) 



for any arbitrary dimensionless function f(x) of the di- 
mensionless Ricci scalar: 
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3oo — n — t ' 

& 

922 = -r 2 , 



9n = -1- 



333 



-r 2 sin 2 9. 



(5) 



for a Newtonian gravitational potential <$>. 

For a particle on circular motion about the mass M 
in the weak field limit, the lowest order of the theory is 
obtained when the left-hand side of equation Q is of or- 
der v 2 /c 2 and when the right-hand side is of order <fi/c 2 . 
Both are just orders C(l/c 2 ) of the theory, or simply 
0(2). As such, when lower or higher order corrections of 
the theory are introduced we will use the notation 0(n) 

re- 



for n = 0, 1, 2, . . . meaning 0(0), ©(c^ 1 ), < 
spectively. 

Having in mind further astrophysical applications (of 
motion of material particles and bending of light - 
massless particles), we expand the metric g^ u about a 
flat Minkowski metric r]^ u :— diag(l, — 1, — 1, — 1) up to 
the second order in time and radial position in such a 
way that 



floo 
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(0) 

=3oo 

J0) 

-3n 

(0) 
=322 



=322 sin 6, 



(2) 

9oo 

(2) 

9n 



0(4), 
+ 0(4), 



(G) 



where the superscript (p) denotes the order 0(p) at 
which a particular quantity is approximated. From equa- 
tions © it follows that the contravariant metric compo- 
nents are given by 



X :— L M R, 



where R is the standard Ricci scalar and 



T „, - ^r 1 /2/ 1 /2 



(9) 



(10) 



is a length scale with 



GM 



lM ■- 



( GM 
V «o 



1/2 



(11) 



with Im the mass-length scale of the system defined by 
[13], and a = 1.2 x 10~ 10 m/s 2 Milgrom's acceleration 
constant [see e.g. [l4|, and references therein] and ( is 
a coupling constant of order one which has to be cal- 
ibrated through astrophysical observations. This f(\) 
theory was constructed by the inclusion of do as a fun- 
damental physical constant, which has been shown to 
be of astrophysical and cosmological relevance [see e.g. 

& SB MM HE m. 

Following the description of Bernal et al. Q the matter 
action takes its ordinary form: 



Sn 



1 



-gd x, 



(12) 



with £ m the Lagrangian matter density of the sys- 
tem. The null variation of the complete action, i.e. 
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8 (Sf + S m ) — 0, with respect to the metric g^ v yields 
the following field equations: 



f'{x)xn» - -^f{x)g^ 



L 2 M (V M V„ - g^A) f'(x) 
8irGLl 



(13) 



M 



where the Laplacc-Bcltrami operator has been written 
as A := V^V^, the prime denotes the derivative with 
respect to the argument and the energy-momentum ten- 
sor T M „ is defined through the standard relation 6S m = 
— {\/2c)T a p8g a P . Also, in equation (|T5jl . the dimension- 
less Ricci tensor is defined as: 



Xnv ■ L^jR^v , 

where R^ v is the standard Ricci tensor. 
The trace of equations (fl3"]l is 



(14) 



f(x)x-V(x) + 3L 2 M Af'( X ) 



*^T, (15) 



where T := T". 

Bcrnal et at [3| have shown that the function fix) 
must satisfy the following limits: 



f(x) = 




when x ^ 1, 
when x "C 1. 



(16) 



The limit x ^ 1 recovers Einstein's general relativity 
and the condition x *C 1 yields a relativistic version of 
MOND. In this last regime it follows [see Q from the 
trace equation (fT"5)) that: 



/'(x)x-2/(x)«3L2 f A/'(x), (17) 

to all orders of approximation and so the trace (|15p can 
be written as: 



3L 2 M Af(x) 



8ttGL 2 a 



T. 



(18) 



Since we are interested in the field produced by a point 
mass M, then the right-hand side of equations (|T3|) , (fT"5|) 
and f|18[) are null and so the last relation in vacuum can 
be rewritten as: 



fix) = x b - 

In this case, relation (fT9|) is equivalent to 



Af(R) = 0, 



(20) 



(21) 



to all orders of approximation for a power-law function 
of the Ricci scalar 



f(R) = R b 



(22) 



Substitution of the power-law function (|20l) in the null 
variations of the gravitational field's action ([8} in vacuum 
means that 



6St 



- ff d 4 a; = 0, 



(23) 



and so 



I?" 



-5d 4 x = 0. 



(24) 



This equation gives the same field equations as the null 
variation of the action for a standard power-law metric 
f(R) theory (|2"2"|) in vacuum. With this in mind, we can 
follow the standard perturbation analysis for f(R) re- 
stricted by the constraint equation (|21l) needed to yield 
the correct MOND- like limit. Since we are only inter- 
ested in a power-law description of gravity far away from 
general relativity (cf. equation (JT6J) ) , then in what fol- 
lows we use the standard f(R) field equations for vacuum 
as described by Capozziello and Faraoni Q for a power- 
law description of gravity given by equation (|22j) with 
b = 3/2, with the constraint (|2"Tj) . To follow their nota- 
tion, we write the field equations (| 13[) in vacuum as 



(25) 



f'(R)R tiV --f(R)g fl v + H^ = 0, 



where the fourth-order terms are grouped into the fol- 
lowing quantity: 



:= - (V„V„ - g^A) f'{R) . (26) 
The trace of equation (|23|) is thus given by 



Af'(x) = 0- (19) f'(R)R-2f(R)+H = 0, (27) 

As shown by Bernal et al. Q, /(%) = x 3 ^ 2 yields the with 
correct MONDian non-relativistic limit. However, for the 
sake of generality we will assume in what follows that the 

function f( X ) is of power-law form H := U^g^ = 3A/'(i?). (28) 
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In what follows the sign convention used in the defini- 
tion of the Riemann tensor becomes a relevant point in 
the subsequent discussion. As discussed in appendix lAl 
the solutions to the differential field equations of any 
f(R) theory of gravity greatly depend on the signature 
chosen for Riemann's tensor, a bifurcation which does 
not appear in Einstein's general relativity. Throughout 
the article we select a particular branch of solutions given 
by the nowadays almost standard definition of Riemann's 
tensor in equation (|A3|) . 

In dealing with some of the cumbersome algebraic ma- 
nipulations that a perturbation to an f(R) theory of 
gravity presents, we have used the Computer Algebra 
System (CAS) Maxima to facilitate the computations. 
The MEXICAS (Metric EXtended-gravity Incorporated 
through a Computer Algebraic System) code (Copyright 
of T. Bernal, S. Mendoza and L.A. Torres and licensed 
with a GNU Public License Version 3) we wrote for this 
is described in appendix [Bl and can be downloaded from: 
|http: / / www.mendozza.org/sergio /mexicas| Further de- 
velopment on the treatment of the field equations by the 
MEXICAS code is described in appendix [Cj 

For the case of a static spherically symmetric space- 
time (p} it follows that 



The fact that — is consistent with the flatness 

of spacetime assumption at the lowest zeroth order of 
perturbation. The expression for the second order com- 
ponent of Ricci's scalar from the metric components © 
is given by 



(2) 

o (2) 



(2) 
%Q,rr 



(2) 



;5oo 



(34) 



The global minus sign that appears on the right-hand side 
of equation p4|) for Ricci's scalar R^ at second pertur- 
bation order differs from that reported by Capozziello 
and Stabile [7], Capozziello, Stabile, and Troisi [8j. As 
mentioned above, and discussed in appendix |A1 this fact 
occurs due to the choice of signs in the definition of Rie- 
mann's tensor. The particular choice used throughout 
the article is the one given by equation (|A3| and so, our 
solutions lie in a different branch as the one reported by 
those authors. 



IV. LOWEST ORDER SOLUTION 



f 5 n (lnV=s) J R,r + g n R 
f"'{ R,tiR,u — Q^vQ R 



9 U r 



11 D 2 
r I i 



and 



n = 3/" [(g n r + g u (in V=g) >r ) R,r + 9 n R,rr 



V'"g ll R,l 



We note that since 



-g = r 2 sin 6 jl 



(2) (2) 



1 1/2 

+ 0(4)} , (31) 



(29) 



(30) 



then, by using the fact that In (y/—g) _ = (y/—g) r I \/~ 9i 
it follows that 



Let us now calculate the order of the trace equa- 
tion (|2"7|) using relations (1221 and (1551) . On the one hand, 
the lowest order of the first two terms on the left-hand 
side of the trace equation is 0(26). On the other hand, 
direct inspection of the right-hand side of equation (|3"01 
results in the fact that the lowest order of T~L is 0(2b— 2). 
Indeed, the last term of the right-hand side of this equa- 
tion is oc 2 and so, to lowest order of pertur- 
bation of relations (0 and (|33|) , this means that T~L con- 
tains terms of the form R^ b ~ 3 R^r and so, H is of order 
0(2b — 2). This analysis indicates that to lowest order 
the trace equation to consider is 



7^(26-2) =3A/ /(26-2) (J?)=0 



(35) 



This result is consistent with relation ([2Tj) to lowest order 
of approximation and is in perfect agreement with the 
perturbative study performed by Bernal et al. [3j. Note 
also that this is the only independent equation at this 
order. 

Direct substitution of equations (|22|) and (|33| into the 
last equation leads to 



In 



f =9). 



2 1 

r + 2 



(2) 

5oo, 



(2) 

9u T 



+ 0(4) . (32) 



Since Ricci's scalar depends on the metric components 
and their derivatives up to the second order with respect 
to the coordinates, it follows it can only have a non-null 
second and higher perturbation orders, i.e. 



R = R (2) + R (4) + 0(6) 



(33) 



rt 2b - 2) = 3&(&-i)fl< 2 > 6 -y i(0) 



(In- 



>(o) 



RV 



+R^ 2 l\ + 36(6 - 1)(6 - 2)i?( 2 > b -y i (°)i?( r 2 ) 2 - 0. 

(36) 

Substitution of expressions (0 and ([32]) in the previous 
equation leads to the following differential equation for 
Ricci's scalar at order 0(2): 
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where 



r ' 



+ (6 - 2)i?< 2)2 = . 



which can be written in a more suitable form as 



lni?| r 2 ) 



(6-2) 



lni?( 2 ) 



The solution of the previous equation is: 



(37) 



(38) 



= ~ (V„ V, - g^A) f'W (R) . (43) 

The complete R.ffi from equation (|29|) is written in ap- 
pendix o 

Now, from equation (|2"Tj) it follows that the Laplace- 
Bcltrami operator applied to f'(R) must be zero at all 
perturbation orders. In particular A/'( 26 ) = 0. With 
this condition, the field equations (|42|) simplify greatly 
and can be written as 



R^(r) = 



ib- 1,(4 

r 



B 



1/(6-1) 



(39) 



where A and B are constants of integration. 

Far away from the central mass, spacetime is flat and 
so Ricci's scalar must vanish at large distances from the 
origin. This means that the constant B — and so 



R^(r) = 



(6-1 



.4 



1/(6-1) 



(40) 



As explained by Bernal et al. [3j , the case 6 = 3/2 yields a 
MOND-likc weak field limit and so, substituting 6 = 3/2 
in relation (l4Tfl) yields: 



*< 2 >(r) = |, 



(41) 



where R := A 2 / 4. This is exactly the same result as the 
one obtained by Bernal et al. 0]. As these authors have 
shown, this result yields a MONDian-like behaviour for 
the gravitational field in the limit r Im 3> r g- For this 
particular case, the lowest order of approximation of the 
theory is 0(1), which has a higher relevance as compared 
to the order 0(2) of standard general relativity for which 
6=1. Using very general arguments, the authors also 
showed that the constant R oc t & /Im is proportional to 
the square root of the mass of the central object. In 
order to calculate R from perturbation analysis we need 
to find the expressions for the metric at order 0(2) of 
approximation. 



-6(6 - 1) ji?( 2 > 6 - 2 [j$ - rjc°> j#> - rj£> jz< r 2 > 

+ (6- 2)R^ b - 3 R^ \r$ - r^ 0) i? ( ? 



(44) 



-6(6- 1)(6- 2) 



2R 



(2)6-3 d(2) d(4) 



+(b-3)R (2)b - 4 R {4) R < i ) R { v ) 



Direct substitution of the following Christofiel symbols 



r 1(0) -o r 1(Z > — -: n u w fl w , i-> 

1 00 1 oo — 2 5 300,r ' K^ ) 

and relations © and (J7J in the 00 component of equa- 
tion (HU) leads to 



,1(2) 



,11(0) J 2 ) 



6i?( 2 ) b - 1 J R 2 ) - ii?( 2 ) 6 + h(b - l)gg] r R^ b -^ = 0. 

(46) 

If we now substitute 6 = 3/2, expression (|4Tj) and the 
value of Ricci's tensor at 0(2) of approximation: 



R 



(2) , o (2) 
(2) _ r 500,rr + / .9oO,r 

" 27 1 



(47) 



into equation (|4l)|). we obtain the following differential 
equation for : 



V. f{x) = X 3/2 METRIC COMPONENTS 

Let us now solve the field equations at the next order 
0(26) of approximation. At this order we expect the 

metric components g$ , g$ and Ricci's scalar R^ to 
play a relevant role in the description of the gravitational 
field. In fact, the field equations at this order are given 
by 

bR^R$ \R^) + H$> = 0, (42) 



^g rr . + 3r 5 g r + — = 0, (48) 

and so 

»»M = -f"'(£) + &' < 49 > 

where k\ and r* are constants of integration. By sub- 
stitution of this result in equation (|3"4"]l and using equa- 
tion (|4"Tj) we get the following differential equation for 

(2) 

9xi- 
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(2) 
r 9u,r 



(2) 

•9u 



r2 + 3 ' 



with solution: 



fc 2 _ i? 
r 3 

where fc 2 is a constant of integration. 



(50) 



(51) 



VI. METRIC COEFFICIENTS FROM 
ASTRONOMICAL OBSERVATIONS 

In this section we derive the constraints that the well 
established astrophysical phenomenology of asymptoti- 
cally flat galactic rotation curves satisfying the Tully- 
Fisher relation, and the cumulative gravitational lensing 
observations for elliptical and spiral galaxies and galaxy 
groups, imply for the metric coefficients for static, spheri- 
cally symmetric spacctimcs for any metric theory of grav- 
ity where dark matter is not required. 

To begin with, let us take the radial component (J3J 
of the geodesic equations (|2J) in the weakest limit of the 
theory. In this limit, the rotation curve for test particles 
bound to a circular orbit about a mass M with circular 
velocity v(r) given by equation @ is 



v 2 (r) 



1 



,ii. 



- "<7 300, r- 

err L 



(52) 



Except for the inner regions of spiral galaxies, v(r) can 
be well approximated by a constant which scales with 
the fourth root of the total baryonic mass Mb of the 
spiral galaxy in question, as described by the Tully-Fisher 
empirical relation [see e.g. 



v = {GM h a Q f'\ 



(53) 



In fact, it is from numerous observations of galactic rota- 
tion curves and total baryonic mass estimates, that the 
constant has been calibrated [see e.g. [13, and refer- 
ences therein]. 

We now substitute equations (JBJ) and (O to order 0(2) 
of approximation and the relation (I53p in equation (152 

(2) 

to obtain the following differential equation for (?q : 



(2) 
%0.r 



having as solution 



2(GM b a ) 1 /2 



(54) 



In 



(55) 



where is a scale radius which, from the point of view 
only of the flat rotation curves of galaxies and the Tully- 
Fisher relation, remains arbitrary. We therefore see that 
a necessary and sufficient condition in any metric rel- 
ativistic theory of gravity, where all observational con- 
straints of galactic rotation curves are satisfied without 

(2) 

invoking dark matter, is that g^ must satisfy the previ- 
ous empirically derived relation. 

Comparing the theoretical metric coefficient 5q given 
by (|49|) (obtained from perturbation theory for fix) — 
X 3 / 2 ) and the empirical one (|55|) (obtained from the phe- 
nomenology of flat rotation curves and the Tully-Fisher 
relation), give the following values for the integration 
constants needed in equation 



fci = 0, R = 6r e /l 



M- 



In this case, the gravitational potential 
tion (O takes the form: 



(56) 
from equa- 



V In 



(GM b a ) 1/2 In 



J* J \ r * 
which yields a radial MONDian acceleration: 



(57) 



\V<f>\ = 



(GM h a ) 1/2 



(58) 



Thus, in the v/c -C 1 limit, the fix) = X 3 ^ 2 presented 
is seen to agree to the observed phenomenology of the 
observed galactic rotation curves in the absence of dark 
matter, as already shown by Bernal et al. 

The gn metric coefficient will be obtained from gravi- 
tational lensing phenomenology. We begin from the gen- 
eral deviation angle equation [see e.g. [33, [43| 



p = 2 



[-goo(r)flii(Q] 1/2 dr 
n r [(r/r l ) 2 3oo(r l ) - 5ooM] 1/2 



(59) 



where r"j is the impact parameter of the problem. 

Over the last few years it has become clear that the 
complete phenomenology of gravitational lensing, at the 
level of extensive massive elliptical galaxies [see e.g. 

groups [see e.g. [32|, clusters of galaxies 
[se e e.g. I25l l33j and more recently spiral galaxies [see e.g. 
[Tl],[4l[ can be accurately modelled using total matter dis- 
tributions having isothermal profiles, when treating the 
problem from the point of view of Einstein's general rel- 
ativity. All these observations show that the dark mat- 
ter halos needed to explain gravitational lensing under 
Einstein's general relativity obey the same Tully-Fisher 
scaling with total baryonic mass as the ones needed to 
explain the observed rotation curves of spiral galaxies. 
This means that for a given total baryonic mass, spiral 
and elliptical galaxies, as well as clusters and groups of 
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galaxies require dark matter halos having the same phys- 
ical properties to explain the observations; from kinemat- 
ics of rotation curves in the former case to gravitational 
lensing in the latter one [IH|4l[. Under Einstein's general 
relativity the majority of these isothermal matter distri- 
bution, particularly at large radii, must be composed of 
a hypothetical dark matter. 

For a static spherically symmetric total matter distri- 
bution Mt, since assuming the validity of Einstein's gen- 
eral relativity Schwarzschild's metric holds, and therefore 
ffoos = -l/5iis, we get: 



2r 2GMtM 
ffoos — -I — 1 5 — 1 — 2 

r err 



(!)"■ « 



The subscript S identifies the coefficients of the 



Schwarzschild metric, and Mr(r) 



r /G refers to the 



hypothetical isothermal total matter distribution [cf. |4[ 
needed to explain the observed lensing, when assuming 
general relativity. From this it follows that the dark 
matter hypothesis provides a self-consistent interpreta- 
tion of observed phenomenology: the same dark matter 
halos, which are required to explain the observed rota- 
tion curves, have been solved for by analysing extensive 
lensing observations. 

From equation (|60|) it follows that for isothermal total 
matter halos under Einstein's general relativity, the met- 
ric coefficient 500s does not depend on the radial coordi- 
nate. We can see this by using the empirical Tully-Fisher 
relation (l53l) between the velocity and the total baryonic 
mass in the last identity above. Thus, the coefficient (|60|) 
can then be taken outside of the integral (|59| of the de- 
viation angle, where for the Schwarzschild metric and 
isothermal total matter halos we now obtain 



[i-2(v/ c y 



i i/2 L 



r[{r/ ri y-l] 



1/2 



7T. (61) 



The above radial integral yields 7r/2 and we obtain the 
observed bending angle as 



[l-2(v/c) 2 } 1 / 2 



[l-2(GM b a ) 1/2 /c 2 ] 1 /2 



(62) 

We see that the well established empirical result of 
lensing observations yielding isothermal total dark mat- 
ter halos under the standard theory is strictly the obser- 
vation of constant bending angles which do not depend 
on the impact parameter, scaling with the observed bary- 
onic total masses as indicated above. 

Now, since (v/c) 2 is of order 0(2) we can write equa- 
tion (|6"2")l as 



the bending angle does not depend on the impact param- 
eter and scales with the square root of the total baryonic 
mass. This last equation gives a clear illustration of the 
link between the dynamics and the spacetime curvature 
effects induced by the presence of an observed baryonic 
mass. 

We can now use the result of equation (j6"3")l to con- 
strain the metric coefficient gn for any metric theory of 
gravity, seeking an accurate description of the observed 
gravitational lensing phenomena without the introduc- 
tion of any hypothetical dark matter. To do this, let 
us return to the general lensing equation (f5§)) , and ask 
that the results obtained under the Schwarzschild met- 
ric with isothermal total matter halos match those under 
any metric theory of gravity, at all impact parameters 
and for any total baryonic masses: 



Q3- 



(- 



dr 



n r[(r/ri) 2 -l 
-9oo(r)gu(r)} 1/2 dr 



1/2 



(64) 



r[{r/ri) 2 go {ri) - g 00 (r)] 



1/2 : 



at 0(2) of approximation from equations (|59[) and (|6T 
Let us rearrange integral (|64|) in such a way that: 



1 



1 

r[(r/ri) 2 
[-goo(r)ffn(r)] 1/2 
r [(r/r l ) 2 .g 00 (r 4 ) - g o(r)] 1/2 



1I/2 



dr = 0. 



(65) 



Since the result must hold for all impact parameters, the 
integrand of the above equation must be equal to zero 
and so, 



u \ 2 
c 



-9oo(r)gu(r) 



(r/r 



- 1 



(r/n) 2 g o(n) - goo(r) 



(66) 

Taking again v/c -C 1 it follows that that the metric 
coefficient gn is given by: 



(r/n) 2 [ffoo(r. t )/ffoo(r)] - 1 
(r/r,) 2 - 1 



From a mathematical point of view, since the contri- 
bution to the integral in the lensing equation ([55)1 is fully 
dominated by the region r»r„ and given the very mild 
radial dependence of the empirical goo term, we can take 
<?oo(fi) ~ 9oo( r ) in the above equation to yield: 



{GM h a o yl 2 r g 

— ~ 2 — = 



A I 



(63) 



The above equation summarises all empirical results of 
gravitational lensing at galactic and galaxy group scales: 



2(GM b a 



,1/2 



Im 



(68) 
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Thus, any metric theory of gravity where gu matches 
the above expression in the regime where gravitational 
lenses are observed will accurately reproduce all the ob- 
served lensing phenomenology, with the total baryonic 
mass of the object in question (galaxies or group of galax- 
ies), and no hypothetical dark matter is assumed to exist. 
Equations (|55p and (|68[) give empirical mathematical re- 
lations for the metric coefficients at perturbation order 
0(2) which reproduce all observed rotation velocity and 
gravitational lensing phenomenology, without the inclu- 
sion of any dark matter component. 

Notice that the mass dependence of the second term 
on the right-hand side in expression (l68l) for the metric 
coefficient gu is the same as the factor in expression (|55|) 
for goo- This last was obtained for a rigorously flat rota- 
tion curve in accordance with the Tully-Fisher relation. 
This shows that the ratio t^/Im of the two important 
characteristic lengths of the extended metric theory of 
gravity proposed by Bernal et al. 3 is the determinant 
dimensionless measure of deviations from flat spacetime 
at galactic scales, exactly as exp ected from the dimen- 
sional analysis in Hernandez [17} 

The metric coefficient gu in equation (|68[) can be di- 
rectly compared to the results for the fix) — X 3 ^ 2 metric 
theory of Bernal et al. [|| obtained in equation (f5Tj) with 
the inclusion of the results of equation ([55| . This means 
that the choice of the integration constant 

fc 2 = 0, (69) 

makes these expressions for the metric component gu 
identical. 

Use of the mathematical approximation A x s» 1+x In A 
to write the following expressions for the full empirical 
metric coefficients gives: 

goo « 1 + (2r g // M ) In (r*/r) « (r*/r) 2rg/ ^ , (70) 
gu*-l-(2r s /l M )*e 2r */ 1 **. (71) 

We note that all the approximations used in this 
section introduce an error several orders of magnitude 
smaller than the intrinsic observational uncertainties in 
the empirical relations used. Therefore, all of the ex- 
pressions given can be considered as strictly equivalent 
in regards to the accurate modelling of astrophysical ro- 
tation curves and gravitational lensing data. 

VII. DISCUSSION 

Through the use of the weak field limit of the met- 
ric f(xL= X 3 ^ 2 theory of gravity constructed by Bernal 
et al. [3(, we have shown that it is possible to explain 
both the dynamics of massive particles and the deflec- 
tion of light by observed astronomical systems such as 
elliptical galaxies, spiral galaxies and groups of galaxies. 
Recently, the same metric theory of gravity was shown 



to be coherent also with the expansion dynamics of the 
observed universe [§:, 26] . This is an expected result from 
a theory of gravity constructed through astronomical ob- 
servations: it must be coherent at all scales. The regime 
of Einstein's general relativity is by no means violated, 
since the applications developed in this article (r > Im) 
lie far away from the mass and length scales associated 
to the ones of Einstein's general relativity (r <C Im) [see 
e.g. EH. 

The results of this article were constructed using a 
static spherically symmetric metric with the time and 
radial components perturbed up to order 0(2) of approx- 
imation. This work generalises the one of Bernal et al. Q 
in which the radial metric component was assumed up to 
order 0(0) only and so, information on the choice of sig- 
nature of the Riemann tensor was lost (see appendix \^ . 
Such information is very important while working with 
fourth order metric theories of gravity. 

Wc mention again the tremendous importance of a cor- 
rect choice for the signature of the Riemann tensor as 
described in appendix [S] The choice (IA3|) . and only that 
choice, used in this article yields results in agreement 
with astronomical observations. In other words, astro- 
nomical observations fix the correct (and unique) choice 
of signature for Riemann's tensor. This is an important 
result, since otherwise solutions from the other branch 
appear which are not in accordance with astronomical 
observations. 

Table U summarises our main results. It is important to 
note that the empirical values of the metric components 

(2) (2) 

3 q and g\^ do not depend on any gravitational theory 
and as such, they represent functions that any successful 
theory of gravity (such as the one used in this article) 
needs to match. 

An important fact arises from the usage of the fix) 
metric theory of gravity and not the f{R) formalism. Al- 
though closely related to each other, the correct dimen- 
sional approach fix) introduces mass and length scales 
that, as shown by Bernal et al. Q, need to be incorpo- 
rated into the gravitational field action. Although the 
field equations in vacuum for both fix) and f(R) under 
a power-law representation yield the same field equations 
(since the mass M generating the gravitational field is a 
constant), /(-R) gravity is not capable of reproducing the 
crucial lensing observations as it lacks a crucial constraint 
equation. The gravitational theory fix) = x 3 ^ 2 is able to 
do so since under this approach the correct limit where 
MONDian-like effects are expected yield the constraint 
equation (|19l) or ([2~Tj) . Notice however that both /(i?) 
and f{x) with the appropriate choice of Riemann's ten- 
sor (|A3|) are able to reproduce the flat rotation curves of 
galaxies and the Tully-Fisher relation. 

In an effort to generalise and look for a fundamental 
basis to an fix) theory of gravity, Carranza, Mendoza, 
and Torres @ and Mendoza [26| have shown that these 
metric theories are equivalent to the the F(R,T) con- 
struction of Harko et al. [l6[. These authors have also 
shown that the particular theory fix) — X 3 ^ 2 is in excel- 
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(2) 

9ao 


(2) 

Sir 


Observations 


l M V r * ) 

(Tully-Fisher) 


-Hit 
(lensing) 


Theory 

fix) = x 3/2 


R = 6r g /l M 


) + * 
fci =0 


T 3 " "r r 3 

R = 6r s /lM fe = 



TABLE I. The table shows the results obtained for the met- 
ric components and for a static spherical symmetric 
spacetime in scales of galaxies and galaxy groups obtained 
empirically from astronomical observations of these systems 
and the ones predicted by the metric f(x) = X 3//2 theory of 
gravity of Bernal et al. A good metric theory of gravity 
must be such that it converges to the inferred values presented 
in the table. The theory /(%) = x 3 ^ 2 i s m perfect agreement 
with the observed metric components. The dimensionless ra- 
tio formed by the quotient of the gravitational radius r g to the 
mass- length scale Im (see equation (|11|) ]) is the determinant 
dimensionless quantity of the problem. Since the metric com- 
ponents determine the "gravitational potential" of the sys- 
tem, the length r* is undetermined. However, since the nat- 
ural length scale of the system is Im one can always asummc 
r* — Im, which also ensures no sign change in the potential 
in equation (|57[) over the domain of applicability r > Im- 



lent agreement with cosmological observations of SNIa. 



An f(x) theory of gravity satisfying the limits of equa- 
tion fp~6|) implies that gravity is no longer scale-invariant. 
In fact, precise gravity tests have been performed only at 
strong regimes of Einstein's gravity, where x » 1, and 
so the involved accelerations of test particles are such 
that a >• a [see e.g. 45]. In exactly the opposite regime, 
where x <§; 1, where the involved accelerations of test 
particles are such that a < ao, gravity differs from Ein- 
stein's general relativity. The traditional approach of 
assuming Einstein's general relativity to be valid at all 
scales means that unknown dark entities are needed to 
explain various astrophysical observations. This article 
heavily reinforces many others 0, 0, i, [jj-lil [H, E3] 
that show how astrophysical and cosmological observa- 
tions can be accounted for without assuming the exis- 
tence of dark entities and extending gravity so as to be 
non scale-invariant. 
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Appendix A: Comments about the sign convention 
in Riemann's tensor 



In the study of the gravitational field equations, the 
link between the curvature of spacetime and the matter 
content is a key fact. All the information regarding the 
curvature of spacetime is contained in the Riemann cur- 
vature tensor R a p v g, which is a function of the first and 
second derivatives of the metric. From a purely math- 
ematical point of view, the Riemann tensor can be ob- 
tained from the commutator of covariant derivatives [l(| : 



{V^V„}Ve = R?V° 



(Al) 



for any vector field V a . From a geometrodynamical 
point of view, the curvature tensor is constructed through 
the change AA^ in a vector after being displaced 
about any infinitesimal closed contour [241 ]: AA^ = 
§T ) ^ lu A\Ax v . By the use of Stokes' theorem it then fol- 
lows that for a sufficiently small closed contour: 



AA„ 



1 



R\ ve A x Af 



(A2) 



where A/" 8 represents the infinitesimal area enclosed by 
the contour of the line integral. In this respect, it fol- 
lows that the Riemann tensor measures the curvature of 
spacetime [cf . |24| . 

In equations (|A1|) and (IA2[) . the Riemann tensor has 
been defined as: 



F A 



(A3) 



If Riemann's tensor is defined by equation (|A3[) . then 
Ricci's tensor is R va := g^^Rp^ua and Ricci's scalar is 
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R a a . Since these are the most used definitions in rela- 
tivity theory nowadays, we will refer to these quantities 
as "standard" . 

However, there is another way in which Riemann's ten- 
sor (and Ricci's tensor) can be defined, usually adopted 
by mathematicians and by Computer Algebr a System s 
(CAS) such as Maxima ( http: / / maxi ma7sourceforge.net[ ) . 
In these cases, the svntaxis is such that [see e.g.l42j| 



= r- 3 — r@ 



p/3 y\ 



(A4) 



If Riemann's tensor is defined by equation (|A4|) . then 
Ricci's tensor is R va :— g^^Rp^va and Ricci's scalar is 
R a a . Although this choice of signs for the Riemann and 
Ricci tensors is not very much in use these days, some 
well-known textbooks use them [see e.g. the Table of 
Sign Conventions at the beginning of reference Hl| . 
The CAS Maxima uses the definition (IA4[) and is such 
that: 



R 



maxima 



R 



■standard j 



(A5) 



in free-index notation. 

As discussed in the Table of Sign Conventions of Mis- 
ner, Thorne, and Wheeler [3l[, general relativity can use 
any of the above definitions (and a few more) simply be- 
cause of the linearity with which Ricci's scalar and Ricci's 
tensor appear in Einstein's field equations. This is how- 
ever not the case in metric f(R) theories of gravity, since 
for example in those theories, the trace of the field equa- 
tions is given by [see e.g. @[: 



f'(R)R-2f(R) + 3Af'(R) 



8ttG 



T. 



(A6) 



To highlight the point, let us substitute the power-law 
function (|22|) in the previous equation to obtain: 



(6 - 2)R b + 36Ai?' 



6-1 



8nG 



T. 



(A7) 



This equation reflects a crucial fact about the choice of 
sign in Riemann's tensor. Due to the presence of the 
derivative term f'(R) — bR b ~ 1 , depending on the sign 
convention of the definition of the Ricci scalar, there ap- 
pears a sign factor (it) 6-1 which is not global to all the 
terms in the equation. This establishes a bifurcation in 
this class of solutions of the theory. Indeed, for a situ- 
ation where f(R) = R a + R b or any more complicated 
function of R, there is not (a priori) any indication of 
which convention in the definition of Riemann's tensor 
should be used to describe a particular physical phenom- 
ena. In this article we show that the convention can 
be settled through the use of astrophysical observations. 



For example, the results presented in this article were 
obtained with the standard definition of Riemann's ten- 
sor in equation (|A3|) . That choice (and only that one) 
can account for both observed dynamics of massive par- 
ticles in spiral galaxies through the Tully-Fisher relation, 
and for the deflection of light observed in gravitational 
lenses. An important aspect to point out is that the case 
f(R) = R of Einstein's general relativity is free from the 
above ambiguity. This is so because it is possible to re- 
define the signature for the energy-momentum tensor to 
recover the same field equations [see e.g. l2ll l3l|. 

We see from this result that previous works by 
Capozziello and Stabile ,7], Capozziello, Stabile, and 
Troisi [|| have selected the convention used by the CAS 
Maxima in order to compute their results. In that re- 
spect, their results lie in another branch of the solutions 
of the field equations. If we would have taken for ex- 
ample, the definition of Riemann's tensor by Maxima, 

(2) 

then the metric coefficients would have been: <7q = 
2R\n 2 (r)/9+A hx(r)+B and gff = -2R\n(r)/9 + D/r + 
(R — A)/2 (where A, B and D are constants). These are 
very different from the ones obtained in equations (|49[) 
and (|5ip and would have never reproduced the astrophys- 
ical observations treated in this article. It is only through 
the correct choice of signs in the definition of Riemann's 
tensor, such as the ones used in the present article and 
represented in equation (| A3|) . that the good agreement 
with the Tully-Fisher relation and lensing observations 
can be correctly obtained. 



Appendix B: Comments about the maxima code 

In this section we give a brief introduction to the 
code we wrote in the Computer Algebra System (CAS) 
Maxima (http://maxima.sourceforge.net) to obtain the 
field equations. Specifically, we work with the module 
ctensor [cf. [42j]. The syntax of such module is that, 
when invoked, it runs an input interface to design the 
form of the covariant metric. 

The Maxima code MEXICAS (Metric EXtended- 
gravity Incorporated through a Computer Algebraic Sys- 
tem) is Copyright of T. Bernal, S. Mendoza and L.A. Tor- 
res, licensed under a GNU Public General License (GPL), 
version 3 (see http://www.gnu.org/licenses) can be ob- 
tained from: http://www.mendozza.org/sergio/mexicas 
(see the section about copyright and usage in that web- 
page). 

For the implementation of the code, we consider a per- 
turbative approach in the parameter e ;= 1/c, such that 
the covariant components of the metric are given by 



goo 
5ii 



+ ^4 



0(4), 
+ 0(4), 



(Bl) 



where the angular components are given by the standard 
expressions for spherical coordinates as shown in equa- 
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tion ([6]). With these equations, it is simple to construct 
the contravariant components of the metric: 



3 00 = 1 - e 2 .^ 4 



0(4), 
+ 0(4). 



(B2) 



With these considerations, the metric is recorded in the 
ctensor module. From this fact, it is simple to invoke all 
the quantities required to construct the field equations, 
either in general relativity or for any extended metric 
theory of gravity. For example, in a descriptive way con- 
cerning the syntaxis of maxima it follows that: 



christof (mcs) 



[1U 5 



and with similar syntaxis for the Riemann tensor, the 
Ricci tensor and the Ricci scalar. 

Due to the fact that the metric has an order parameter 
e, all the tensorial quantities involved in the construction 
of the field equations will gain this dependence. In the 
formalism of the code, it is a crucial fact to extract the 
perturbation order of every metric quantity to construct 
the field equations at the desired perturbation order. For 
example, for a generic quantity q calculated from the ma- 
nipulation of the metric, if we consider that q( n ' repre- 
sents such quantity at order n, we have: 



= Hm q, 
e— >o 



(B4) 



which reproduces the flat spacetime limit. For the second 
order we have 



q^ = lim 

■>o 



q-q 



,(0) -pfStif" 



(B5) 



and consequently the fourth order is obtained by 



Similarly, higher perturbation orders can be obtained by 
the obvious generalisation of the previous relation. 

In equations (|B5|) and (|B6I) . it is implied that the 
first order quantities vanish, as is also the case for the 
Christoffel symbols. This computational procedure gives 
as an output a key result used in the article correspond- 
ing to Ricci's scalar at second perturbation order, given 
by equation (IM)) . 



Appendix C: Extended field equations using Maxima 

By using the Computer Algebra System (CAS) Max- 
ima and the MEXIC AS code (see appendix [Bj , we ob- 
tained the field equations up to the second order. 



The trace (|27|) of the field equations (|42j) to the order 
0(26) of approximation can be simplified with the aid of 
the solutions found at the lowest order of approximation 
in Section HVl to obtain 



?(2) 



(6 - 2)i?( 2 ' 4 - 36(6 - l)i?< 2) 

+ \R { ? (flgr + *Sr)] +2(b-2)R<;?R(?} ^ 
+ 3b(b-l)(b~2)R^ 2 R (4 ^0. 



The components Hffi of the field equations ([251 at 



( B3 ) order 0(26) are given by: 
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+ (b-2)R ( ^ b - 3 R^ 
x ({\n^)fRf+R% 
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- 6(6 - 1)(6 - 2) 



2i*W„ 4) ~ 9 f} 2g^RfRf + g^ 
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Dividing the field equations P2")l by i?( 2 ) b ~ 4 and using 
the trace (|C1|) and the last equation, a reduced expression 
for the field equations is found: 



26-1 
6 



R m + 6^(2)3^(2) _ b[b _ 1)i? (2)2 r R w 



_pl(0)o(4) _ pl(2)o(2 



-6(6-l)(6-2)i? (2) 
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6(6 - 1)(6 - 2)(6 - 3)i? (2) i? (2) i? (4) = , 



(C3) 



which can also be regarded as the traceless component of 
the field equations. 



